Structural Features and Zeolite Stability: A Linearized Equation Approach.

Zeolite stability, in terms of lattice energy, is revisited from a crystal-chemistry point of view. A linearized equation relates the zeolite lattice energy using simple structural data readily available from experiments or modeling. The equation holds for a large range of zeolite energies, up to 3 eV per tetrahedron with respect to quartz, and has been validated internally via two simple machine learning automatic procedures for data fitting/reference partitions and externally using data from recently synthesized zeolites. The approach is certain in locating those recently synthesized zeolites in the energy range of those experimentally known zeolites used in the parametrization of the linearized equation. Hidden intrinsic structural data-energy correlations were found for data sets built from energy-relaxed structures along with energy values computed using the same energy functions employed in the structural relaxation. The asymmetry of the structural features is relevant for an accurate description of the energy.

Nucleation of zeolitic imidazolate frameworks: from molecules to nanoparticles.

We have studied the clusters involved in the initial stages of nucleation of Zeolitic Imidazolate Frameworks, employing a wide range of computational techniques. In the pre-nucleating solution, the prevalent cluster is the ZnIm4 cluster (formed by a zinc cation, Zn2+, and four imidazolate anions, Im-), although clusters such as ZnIm3, Zn2Im7, Zn2Im7, Zn3Im9, Zn3Im10, or Zn4Im12 have energies that are not much higher, so they would also be present in solution at appreciable quantities. All these species, except ZnIm3, have a tetrahedrally coordinated Zn2+ cation. Small ZnxImy clusters are less stable than the ZnIm4 cluster. The first cluster that is found to be more stable than ZnIm4 is the Zn41Im88 cluster, which is a disordered cluster with glassy structure. Bulk-like clusters do not begin to be more stable than glassy clusters until much larger sizes, since the larger cluster we have studied (Zn144Im288) is still less stable than the glassy Zn41Im88 cluster, suggesting that Ostwald's rule (the less stable polymorph crystallizes first) could be fulfilled, not for kinetic, but for thermodynamic reasons. Our results suggest that the first clusters formed in the nucleation process would be glassy clusters, which then undergo transformation to any of the various crystal structures possible, depending on the kinetic routes provided by the synthesis conditions. Our study helps elucidate the way in which the various species present in solution interact, leading to nucleation and crystal growth.

Unravelling the key factors in the chlorine-promoted epoxidation of ethylene over a silver-copper oxide nanocatalyst.

Ethylene oxide is one of the most important raw materials in the chemical industry, with an annual production close to 35 million metric tons. Despite its importance, to date, no metal has been found that can compete with the original silver bulk material catalyst discovered in 1931. Recently, a few copper and copper-silver based nanostructures have demonstrated remarkable selectivity and activity, especially when coupled with an industrial chlorine promoter. The present work evaluates the mechanistic role of chlorine as an active promoter of the selective oxidation of ethylene to ethylene oxide in the presence of a silver-copper oxide hybrid nanocatalyst (AgCuO). Experimental kinetic studies combined with density functional theory (DFT) calculations provide insight into the influence that Ag/CuO-supported chlorine atoms have over the ethylene epoxidation reaction. Remarkably, the typically described indirect route via the formation of an oxametallacycle (OMC) is also accompanied by a direct route. Furthermore, the presence of chlorine seems to facilitate a more favorable adsorption energy for ethylene oxide (EO) than for acetaldehyde (AA), the main reaction by-product. As a result, complete oxidation of EO can be further prevented in the presence of this AgCuO hybrid heteronanostructure.

Screening heteroatom distributions in zeotype materials using an effective Hamiltonian approach: the case of aluminogermanate PKU-9.

We introduce a method to allow the screening of large configurational spaces of heteroatom distributions in zeotype materials. Based on interatomic potential calculations of configurations containing up to two heteroatoms per cell, we parameterize an atomistic effective Hamiltonian to describe the energy of multiple substitutions, with consideration of both short- and long-range interactions. Then, the effective Hamiltonian is used to explore the full configurational space at other compositions, allowing the identification of the most stable structures for further analysis. We illustrate our approach with the aluminogermanate PKU-9, where we show that increasing the aluminium concentration changes the likely siting of Al, in agreement with experiment.

Fitting electron density as a physically sound basis for the development of interatomic potentials of complex alloys.

The development of new interatomic potentials to model metallic systems is a difficult task, due in part to the dependence between the parameters that describe the electron density and the short-range interactions. Parameter search methods are prone to false convergence. To solve this problem, we have developed a methodology for obtaining the electron density parameters independently of the short-range interactions, so that physically sound parameters can be obtained to describe the electron density, after which the short-range parameters can be fitted, thus reducing the complexity of the process and yielding better interatomic potentials. With the new method we can develop self-consistent, accurate force fields, using solely calculations, without the need to fit to experimental data. Density functional theory calculations are used to compute the observables with which the potential is fit. We applied the method to a Ni-based Inconel 625 superalloy (IN625), modelled here as Ni, Cr, Mo and Fe solid solution alloys. The capability of the force fields developed using this new method is validated, by comparing the structural and thermo-elastic properties predicted with the force fields, with the corresponding experimental data, both for single crystals and polycrystalline alloys.

Chemical Engineering of Photoactivity in Heterometallic Titanium-Organic Frameworks by Metal Doping.

We report a new family of titanium-organic frameworks that enlarges the limited number of crystalline, porous materials available for this metal. They are chemically robust and can be prepared as single crystals at multi-gram scale from multiple precursors. Their heterometallic structure enables engineering of their photoactivity by metal doping rather than by linker functionalization. Compared to other methodologies based on the post-synthetic metallation of MOFs, our approach is well-fitted for controlling the positioning of dopants at an atomic level to gain more precise control over the band-gap and electronic properties of the porous solid. Changes in the band-gap are also rationalized with computational modelling and experimentally confirmed by photocatalytic H2 production.

Selective sulfur dioxide adsorption on crystal defect sites on an isoreticular metal organic framework series.

The widespread emissions of toxic gases from fossil fuel combustion represent major welfare risks. Here we report the improvement of the selective sulfur dioxide capture from flue gas emissions of isoreticular nickel pyrazolate metal organic frameworks through the sequential introduction of missing-linker defects and extra-framework barium cations. The results and feasibility of the defect pore engineering carried out are quantified through a combination of dynamic adsorption experiments, X-ray diffraction, electron microscopy and density functional theory calculations. The increased sulfur dioxide adsorption capacities and energies as well as the sulfur dioxide/carbon dioxide partition coefficients values of defective materials compared to original non-defective ones are related to the missing linkers enhanced pore accessibility and to the specificity of sulfur dioxide interactions with crystal defect sites. The selective sulfur dioxide adsorption on defects indicates the potential of fine-tuning the functional properties of metal organic frameworks through the deliberate creation of defects.

Modelling a Linker Mix-and-Match Approach for Controlling the Optical Excitation Gaps and Band Alignment of Zeolitic Imidazolate Frameworks.

Tuning the electronic structure of metal-organic frameworks is the key to extending their functionality to the photocatalytic conversion of absorbed gases. Herein we discuss how the band edge positions in zeolitic imidazolate frameworks (ZIFs) can be tuned by mixing different imidazole-based linkers within the same structure. We present the band alignment for a number of known and hypothetical Zn-based ZIFs with respect to the vacuum level. Structures with a single type of linker exhibit relatively wide band gaps; however, by mixing linkers of a low-lying conduction edge with linkers of a high-lying valence edge, we can predict materials with ideal band positions for visible-light water splitting and CO2 reduction photocatalysis. By introducing copper in the tetrahedral position of the mixed-linker ZIFs, it would be possible to increase both photo-absorption and the electron-hole recombination times.

Critical Role of Dynamic Flexibility in Ge-Containing Zeolites: Impact on Diffusion.

Incorporation of germanium in zeolites is well known to confer static flexibility to their framework, by stabilizing the formation of small rings. In this work, we show that the flexibility associated to Ge atoms in zeolites goes beyond this static effect, manifesting also a clear dynamic nature, in the sense that it leads to enhanced molecular diffusion. Our study combines experimental and theoretical methods providing evidence for this effect, which has not been described previously, as well as a rationalization for it, based on atomistic grounds. We have used both pure-silica and silico-germanate ITQ-29 (LTA topology) zeolites as a case study. Based on our simulations, we identify the flexibility associated to the pore breathing-like behavior induced by the Ge atoms, as the key factor leading to the enhanced diffusion observed experimentally in Ge-containing zeolites.

Controlling Thermal Expansion: A Metal-Organic Frameworks Route.

Controlling thermal expansion is an important, not yet resolved, and challenging problem in materials research. A conceptual design is introduced here, for the first time, for the use of metal-organic frameworks (MOFs) as platforms for controlling thermal expansion devices that can operate in the negative, zero, and positive expansion regimes. A detailed computer simulation study, based on molecular dynamics, is presented to support the targeted application. MOF-5 has been selected as model material, along with three molecules of similar size and known differences in terms of the nature of host-guest interactions. It has been shown that adsorbate molecules can control, in a colligative way, the thermal expansion of the solid, so that changing the adsorbate molecules induces the solid to display positive, zero, or negative thermal expansion. We analyze in depth the distortion mechanisms, beyond the ligand metal junction, to cover the ligand distortions, and the energetic and entropic effect on the thermo-structural behavior. We provide an unprecedented atomistic insight on the effect of adsorbates on the thermal expansion of MOFs as a basic tool toward controlling the thermal expansion.

Comparing gas separation performance between all known zeolites and their zeolitic imidazolate framework counterparts.

To find optimal porous materials for adsorption-based separations is a challenging task due to the extremely large number of possible pore topologies and compositions. New porous material classes such as Metal Organic Frameworks (MOFs) are emerging, and hope to replace traditionally used materials such as zeolites. Computational screening offers relatively fast searching for candidate structures as well as side-by-side comparisons between material families. This work is pioneering at examining the families comprised by the experimentally known zeolites and their respective Zeolitic Imidazolate Framework (ZIF) counterparts in the context of a number of environmental and industrial separations involving carbon dioxide, nitrogen, methane, oxygen, and argon. Additionally, unlike related published work, here all the targeted structures have been previously relaxed through energy minimization. On the first level of characterization, we considered a detailed pore characterization, identifying 24 zeolites as promising candidates for gas separation based on adsorbate sizes. The second level involved interatomic potential-based calculations to assess the adsorption performance of the materials. We found no correlation in the values of heat of adsorption between zeolites and ZIFs sharing the same topology. A number of structures were identified as potential experimental targets for CO2/N2, and CO2/CH4 affinity-based separations.

ε-Keggin-based coordination networks: Synthesis, structure and application toward green synthesis of polyoxometalate@graphene hybrids.

Four coordination networks based on the {ε-PMo(V)(8)Mo(VI)(4)O(40)(OH)(4)Zn(4)} Keggin unit (εZn) have been synthesized under hydrothermal conditions. (TBA)(3){PMo(V)(8)Mo(VI)(4)O(36)(OH)(4)Zn(4)}[C(6)H(4)(COO)(2)](2) (ε(isop)(2)) is a 2D material with monomeric εZn units connected via 1,3 benzenedicarboxylate (isop) linkers and tetrabutylammonium (TBA) counter-cations lying between the planes. In (TPA)(3){PMo(V)(8)Mo(VI)(4)O(37)(OH)(3)Zn(4)}[C(6)H(3)(COO)(3)] (TPA[ε(trim)](∞)), 1D inorganic chains formed by the connection of εZn POMs, via Zn-O bonds, are linked via 1,3,5 benzenetricarboxylate (trim) ligands into a 2D compound with tetrapropylammonium (TPA) cations as counter-cations. (TBA){PMo(V)(8)Mo(VI)(4)O(40)Zn(4)}(C(7)H(4)N(2))(2)(C(7)H(5)N(2))(2)·12H(2)O (ε(bim)(4)) is a molecular material with monomeric εZn POMs bound to terminal benzimidazole (bim) ligands. Finally, (TBA)(C(10)H(10)N(4))(2)(HPO(3)){PMo(V)(8)Mo(VI)(4)O(40)Zn(4)}(2)(C(10)H(9)N(4))(3)(C(10)H(8)N(4)) (ε(2)(pazo)(4)) is a 1D compound with dimeric (εZn)(2) POMs connected by HPO(3)(2-) ions and terminal para-azobipyridine (pazo) ligands. In this compound an unusual bond cleavage of the central N[double bond, length as m-dash]N bond of the pazo ligand is observed. We report also a green chemistry-type one-step synthesis method carried out in water at room temperature using ε(2)(pazo)(4) and ε(isop)(2) as reducing agent of graphite oxide (GO) to obtain graphene (G). The POM@G hybrids were characterized by X-ray photoelectron spectroscopy, Raman spectroscopy, powder X-ray diffraction, energy dispersive X-ray analysis, infrared spectroscopy, scanning electron microscopy, transmission electron microscopy and cyclic voltammetry.

Mg/Ca partitioning between aqueous solution and aragonite mineral: a molecular dynamics study.

We have calculated the concentrations of Mg in the bulk and surfaces of aragonite CaCO(3) in equilibrium with aqueous solution, based on molecular dynamics simulations and grand-canonical statistical mechanics. Mg is incorporated in the surfaces, in particular in the (001) terraces, rather than in the bulk of aragonite particles. However, the total Mg content in the bulk and surface of aragonite particles was found to be too small to account for the measured Mg/Ca ratios in corals. We therefore argue that most Mg in corals is either highly metastable in the aragonite lattice, or is located outside the aragonite phase of the coral skeleton, and we discuss the implications of this finding for Mg/Ca paleothermometry.

Zeolitic polyoxometalates metal organic frameworks (Z-POMOF) with imidazole ligands and epsilon-Keggin ions as building blocks; computational evaluation of hypothetical polymorphs and a synthesis approach.

We investigate here a new family of zeolitic Metal Organic Frameworks (MOFs) based on imidazole (im) as the ligand and epsilon-type Keggin PolyOxoMetalates (POMs) as building units. The POM used in this study is the epsilon-{PMo(12)O(40)} Keggin isomer capped by four Zn(ii) ions (noted epsilon-Zn) in tetrahedral coordination. We describe here our methods to first construct and then evaluate the stability of hypothetical 3-D POMOFs possessing a tetrahedral network, typified by dense silica polymorphs and zeotypes and referred here to as Z-POMOFs. We use the analogy between the connectivity of silicon ion in dense minerals or zeolites and the epsilon-Zn, using imidazolate ligands to mimic the role of oxygen atoms in zeolites. Handling the epsilon-Keggin and imidazole as the constitutive building-blocks, a selection of 40 polymorphs were constructed and their relative stabilities computed. Among these Z-POMOFs, the cristobalite-like and zni-structure were identified as the most stable candidates. In parallel, we have attempted to synthesize Z-POMOF structures with epsilon-Zn POMs, synthesized in situ under hydrothermal conditions, and imidazole ligands. We present our first experimental result, the extended material [NBu(4)][PMo(V)(8)Mo(VI)(4)O(37)(OH)(3)Zn(4)(im)(Him)], named epsilon(im)(2). The structure of the hybrid framework is built by the connection of dimerized epsilon-Zn POMs to imidazole ligands in two directions. The obtaining of the first POMOF based on imidazole ligand is an encouraging step towards the synthesis of a new family of POMOFs.

Zeolitic polyoxometalate-based metal-organic frameworks (Z-POMOFs): computational evaluation of hypothetical polymorphs and the successful targeted synthesis of the redox-active Z-POMOF1.

The targeted design and simulation of a new family of zeolitic metal-organic frameworks (MOFs) based on benzenedicarboxylate (BDC) as the ligand and epsilon-type Keggin polyoxometalates (POMs) as building units, named here Z-POMOFs, have been performed. A key feature is the use of the analogy between the connectivity of silicon in dense minerals and zeolites with that of the epsilon-type Keggin POMs capped with Zn(II) ions. Handling the epsilon-Keggin as a building block, a selection of 21 zeotype structures, together with a series of dense minerals were constructed and their relative stabilities computed. Among these Z-POMOFs, the cristobalite-like structure was predicted to be the most stable structure. This prediction has been experimentally validated by the targeted synthesis of the first experimental Z-POMOF structure, which was strikingly found to possess the cristobalite topology, with three interpenetrated networks. Crystals of [NBu(4)](3)[PMo(V)(8)Mo(VI)(4)O(36)(OH)(4)Zn(4)(BDC)(2)].2H(2)O (Z-POMOF1) have been isolated under hydrothermal conditions from the reduction of ammonium heptamolybdate in the presence of phosphorous acid and Zn(II) ions. Tetrabutylammonium cations play the role of counterions and space-filling agents in this tridimensional interpenetrated framework. Moreover, the electrochemistry of the epsilon-Keggin POM is maintained and can be exploited in the insoluble Z-POMOF1 framework, as demonstrated by the electrocatalytic reduction of bromate.

Interplay of water, extra-framework cations and framework atoms in the structure of low-silica zeolites: the case of the natural zeolite Goosecreekite as studied by computer simulation.

Computational methods are described that model accurately the structure of hydrated Ca-bearing zeolites. Using Goosecreekite as a model system we probe the influence of framework ordering, cation siting and hydration of pores on the structure and its stability. We develop a methodology which allows the location of Al within the framework to be determined together with the position of extra-framework cations, in a stepwise fashion, progressing from an anhydrous model, via a dielectric continuum model, to finally, a fully atomistic model of the water in the intrazeolite pore space. Our methods reveal the complex interplay of short- and long-range interactions on the optimal structure of such materials.